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I . INTRODUCTION 


Tliere  is  a continuing  need  for  modeling  the  structural  damage 
inflicted  by  mine  blast  on  the  undersides  of  tanks  and  other  armored 
vehicles,  for  both  vulnerability  assessments  and  to  rationally  provide 
hardening  for  such  vehicles.  This  has  been  a subject  of  long-standing 
interest  at  the  Ballistic  Reseai*ch  Laboratory  and  at  other  agencies  at 
Aberdeen  Proving  Ground.  A considerable  number  of  tests  have  been  per- 
formed and  much  field  damage  data  gathered  over  the  years^'^®.  In  the 
recent  past,  an  analytical  approach  has  been  reported^ 


^Hanna,  J.  V. j "An  Effectiveness  Evaluation  of  Several  Types  of  Antitank 
Mines,  " BHL-MR-616,  US  Army  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  MD  (June  1962).  (AD  #377342) 

"^Bailey,  R.  A.,  Bom,  D.,  atvi  Sultanoff,  M.,  "Analysis  of  the  Performance 
of  the  Mock-Up  Booster  Assembly  for  the  Multi-Jet,  Shaped  Charge,  Anti- 
Tank  Mine,  " BRL-MR-694 , US  Army  Ballistic  Reseai'ch  Laboratory,  Aberdeen 
Proving  Ground,  MD  (October  1961).  (AD  #377333) 

^Haskell,  D.  F.  and  Reisinger,  M.  J.,  "Armored  Vehicle  Vulnerability 
Analysis  Model  - First  Version,  Introduction,"  BRL-R-1867,  US  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD  (February 
1976).  (AD  #B009638L) 


^ Clark,  E.  L. , "Testing  of  Anti-Armor  Devices,"  BRL-CR-221 , US  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD  (April  1976). 
, (AD  #B003825L) 

^Norman,  U.  M. , "An  Estimate  of  the  Performance  of  the  FRG  Pandora  Antv- 
Tank  Mine"  (U),  BRL-MR-1764,  US  Army  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  MD  (July  1966)  ( SECRET -NOFORN ) . (AD  #376685) 


^Cioffi,  A.  R.  and  Vincent,  A.  R.,  "Preliminary  Estimates  of  the  Vulner- 
ability of  Light  Weight  Armored  Vehicles  to  Attack  by  Antitank  Mines," 
(U),  BRL~TN-1197 , US  Army  Ballistic  Research  Laboratory,  t\herdacn  I'roving 
Ground,  MD  (June  1968)  (CONFIDENTIAL) . (AD  #378697) 

^Norman,  R.  M. , "An  Estimate  of  the  Performance  of  the  XM24  Anti-Tank 
Mine"  (U),  BRL-MR-1774,  US  Army  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  MD  (August  IBs's)  (CONFIDENTIAL) . 

’^Norman,  R,  M.  and  Beiahler,  G.  P.,  "Effectiveness  Studies  of  the  XM34  AT 
Mine"  (U),  BRL-TN-1640,  US  Army  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  MD  (December  1966)  (CONFIDENTIAL) . (AD  #379091) 


^Norman,  R.  M,  and  Beichler,  G.  P, , "An  Estimate  of  the  Performance  of 
Buried  XM34  Antitank  Mines"  (U),  BRL-TN~1641 , US  Army  Ballistic  Research 
Laboratory,  Aberdeen  Pr'uving  Ground,  MD  (December  1966)  (CONFIDENTIAL) . 

^ , "Deformation  in  Flat  Plates  Exposed  to  HE  Mine  Blast," 

AMSAA-TM-74,  US  Army  Mateinel  Systems  Analysis  Agency,  Aberdeen  Proving 
Ground,  MD  (May  1970). 

^^Haskell,  D.  F.,  "Deformation  and  Fracture  of  Tank  Bottom  Hull  Plates 
Sub.jeated  to  Mine  Blast,"  BRL-R-1687 , US  Army  Ballistic  Researc' 
Laboratory,  Aberdeen  Proving  Ground,  MD  (May  1972).  (AD  #901628) 
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More  recently,  preliminary  computations,  using  the  two-dimensional 
Eulerian,  hydrodynamic  computer  code  DORF^^,  of  the  time-history  of  the 
loading  of  a rigid  plate  by  a bare,  spherical  charge  have  been  performed^ 

The  DORF  hydrocode  uses  a centered,  finite  difference  technique  to  solve 
the  unsteady  Euler  equations;  it  has  an  artificial  viscosity  option. 

Those  calculations  were  performed  to  determine  the  applicability  of  the 
DORF  hydrocode  to  this  problem,  and  to  gain  some  insight  as  to  the  effect 
of  the  computational  mesh  on  the  solution.  The  bare  charge  computations 
treat  both  the  air  and  the  e.xplosivc  products  as  polytropic  gases,  with 
the  explosive  being  represented  by  a sphere  of  high-pressure  quiescent 
explosive  products  initially  resting  tangent  to  a reflecting  bottom  boun- 
dary. The  top  boundary  is  also  reflective,  representing  a rigid,  fixed, 
target  plate.  Computational  studies^**  using  the  DORF  code  indicate  that 
the  target  plate  loading  using  this  representation  of  the  explosive  at 
initiation  does  not  materially  differ  from  that  computed  by  using  Taylor's^^ 
solution  for  the  detonation  wave.  The  intent  of  the  present  study  is  to 
use  the  experience  gained  from  the  previous  work’^>^'^  to  attempt  to  com- 
pute the  target  plate  loading  and  deformation  due  to  the  blast  from  a land 
mine  buried  in  soil,  thereby  also  testing  the  ability  of  the  DORF  code  to 
model  a soil  cover. 

A recent  study' ^ at  the  Southwest  Research  Institute,  hereafter  re- 
ferred to  as  SwRI,  involves  the  collation  of  existing  data  for  scaled 
experiments  involving  the  encounter  of  buried  land  mines  with  various 
targets,  and  also  reports  new  data  for  a series  of  small  scale  experi- 
ments, involving  flat  .plates,  performed  by  .SwRI.  A correlation  function, 
based  on  the  use  of  dimensional  analysis,  for  the  total  impulse  delivered 
by  the  land  mine  and  overburden  on  a target  is  presented.  For  the  special 
case  wherein  the  target  is  an  initially  flat  plate,  the  report  also  pre- 
sents a function  which  relates  the  strain  energy  stored  in  the  plate  to 
the  total  energy  delivered  to  the  plate  by  the  mine  blast  to  compute  the 
maximum  residual  deformation  of  the  plate.  This  presents  an  opportunity 
to  compare  the  integrated,  experimental,  total  impulse,  as  computed  using 
the  SwRI  correlation  function,  with  that  computed  using  the  DORF  hydrocode. 


^‘^Johnson,  Wallace  E.,  "Code  Correlation  Study,:  AFWL~TE-70-l 44 , Air 
Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  NM  (April  1971). 

^^Miller,  James  E.,  "Preliminary  Study  of  Target  Load  Prediction  by  Use 
of  a Hydrodynamic  Computer  Code,"  BHL  MR-2472,  US  Army  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  MD  (April  1975).  (AD  B00.1829I,) 

Allison,  W.  D.,  BRL,  private  communication  of  as  yet  unpublished  data. 

^^Taylor,  Geoffrey  I.,  Sir,  The  Dynamics  of  the  Combustion  Products 
behind  Plane  and  Spherical  Detonation  Fronts  in  Explosives,  Proceedings 
of  the  Royal  Society  of  London,  Series  A,  Vol.  200,  22  February  1950. 

^^Westine,  Peter  S.  and  Hokanson,  James  C.,  "Procedures  to  Predict  Plate 
Defc'rmations  from  land  Mine  Explosions"  (U) , TACOM  Technical  Report 
No.  12049,  US  Army  Tank  Automotive  Cormand,  Warren,  Ml  (August  1975) 
(CONFIDENTIAL) . 


10 


Maximum  residual  deformation  data  for  the  SwRI  test  series  using  flat 
plates  are  also  reported^®,  as  is  a function  for  the  general  shape  of  the 
deformed  plate.  The  loading  generated  by  the  DORF  hydrocode  is  used  as 
input  to  the  REPSIL^^  structural  response  code  to  compute  the  target  plate 
response.  The  computed  response  is  then  compared  directly  to  the  maximum 
residual  deformation  and  the  general  deformed  shape  which  are  predicted  by 
the  SwRI  correlation  function. 


II.  SIMULATED  SwRI  TEST 

The  SwRI  data  which  show  the  least  amount  of  scatter  on  the  impulse 
correlation  curve  are  the  bulge  plate  tests.  The  SwRI  bulge  plate  test^^ 
selected  for  computer  simulation,  test  number  36  in  the  te.st  series  for 
flat,  square  plates  with  a charge  buried  directly  below  the  plate  center, 
involves  a square  armor  plate,  25.4  cm  on  a side  and  0.678  cm  thick, 
centered  above  a 0.265  kg  charge  of  C-4.  The  total  stand-off  distance, 
the  distance  from  the  bottom  of  the  target  plate  to  the  center  of  the 
charge,  is  16.2  cm,  of  which  8.9  cm  is  the  air  gap  between  the  overburden 
and  the  target  plate.  In  the  experiment,  the  target  plate  is  restrained 
by  a heavy  collar;  however,  the  collar  is  not  massive  enough  to  prevent 
the  plate  and  collar  from  being  thrown  several  feet  in  the  air  by  the 
exploding  land  mine. 

The  high-explosive  C-4  charge  is  of  a pancake  shape,  with  an  assumed 
thickness-to-diameter  ratio  of  5/13.  The  thickness-to-diameter  ratio  for 
the  charge  used  in  this  test  was  not  reported^^.  There  were  references  in 
the  text  to  this  5/13  ratio,  and  a photograph  of  a pancake  charge  to  be 
used  for  a bulge  plate  test  indicated  such  ratio  would  be  a reasonable 
assumption.  Assuming  a normal  density  of  1.592  g/cm^,  the  radius  of  the 
charge  is  4.1  cm,  the  thickness  3.2  cm.  The  distance  from  the  soil  sur- 
face to  the  top  of  the  charge  is  5.7  cm. 

The  geometry  of  the  problem  allows  the  use  of  a cylindrical  coordinate 
system  in  the  hydrocode  calculation.  The  flow  field  prior  to  the  detona- 
tion of  the  charge  is  shown  in  Figure  la.  The  left  boundary  of  the  flow 
field  is  the  axis  of  symmetry.  The  right  and  bottom  boundaries  are  trans- 
missive boundaries,  far  enough  removed  from  the  charge  and  the  top  boundary 
so  that  artificial  reflections  from  these  boundaries  do  not  have  time  to 
invalidate  the  near-field  solution  at  the  top  boundary.  The  computation 
was  stopped  and  analyzed  frequently,  so  that  any  flow  phenomena  that  did 
not  seem  to  have  a physical  basis  could  be  identified,  and  their  impact 
on  the  solution  evaluated.  The  top  boundary  is  a rigid,  immobile, 
reflecting  boundary,  intended  to  simulate  the  target  plate.  This  repre- 
sentation of  the  target  plate  is  a concession  to  the  DORF  hydrocode, 
replacing  a finite,  movable,  and  deformable  plate-and-col lar  combination 


Santiago,  J.  M.,  Wienieuaki,  H.  L. , and  Huffington,  N.  J,,  Jr,,  "A 
User’s  Manual  for  the  REPSIL  Code,”  BRL  Report  1744,  US  Army  Balliatia 
Research  Laboratory,  Aberdeen  Proiting  Ground,  MD  (October  1974). 

(AD  #A003176) 
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The  flow  field  at  cycle  0,  time  = 0. 
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with  an  infinite,  fixed,  and  rigid  reflecting  boundary.  This  concession 
simplifies  the  computational  flow  field  considerably;  it  is  believed  to 
cause  only  a modest  over-estimate  of  the  forces  on  the  plate,  primarily 
by  eliminating  the  possibility  of  computing  rarefaction  waves  which  would 
exist  under  experimental  conditions. 

The  finite-difference  computational  grid  for  this  use  of  the  DORF 
hydrocode  contains  60  flow  field  cells  in  the  radial  direction,  with  the 
vertical  area  projections  of  the  cells  increasing  in  a geometric  progres- 
sion with  increasing  radius.  There  are  80  cells  in  the  axial  direction; 
the  cell  heights  are  determined  by  a series  of  geometric  progressions 
centered  on  the  high-explosive  charge  (see  Figure  lb).  Because  the  DORF 
hydrocode  uses  a centered-differencing  technique  assuming  equal  cell  vol- 
umes, it  is  assumed  that  the  radial  progression  in  cell  volume  used  for 
this  grid  introduces  less  error  than  a radial  progression  in  differential 
rad  i us . 


III.  EQUATIONS  OF  STATE 

Both  the  air  and  the  explosive  products  are  assumed  to  be  polytropic 
gases^^’^^,  governed  by  the  equation  of  state 

p = p I (y-1)  , (1) 

where  p represents  pressure,  p represents  density,  I represents  specific 
internal  energy,  and  y represents  the  ratio  of  specific  heats.  For  the 
air , 


Yg  = 1.4,  (2.1) 

and,  for  the  explosive  products  * ^ , 

Yp  = 2.71.  (2.2) 


The  soil  surrounding  the  explosive  products  is  governed  by  the  Tillotson 
equation  of  state^®*^”^  as  used  by  Johnson  in  the  DORF*^  code,  where  the 
equation  for  the  pressure  takes  on  three  main  forms.  The  pressure  in  the 
compressed  state,  or  any  cold  state,  is 


a + 


I,.n‘ 


+ 1 


Ip  + Ap  + Bp 


(3.1) 


^^Tiltotson,  J.  H.y  "Metallia  Equations  of  State  for  Hyperveloaity  Impact," 
GA-3216,  General  Atomic,  Division  of  General  Dynamics,  San  Diego,  CA 
(July  li)62). 

L.  H.  and  Johnson,  W.  E.,  "Energy  Partitioning  during  Hyper-velocity 
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where  either 


n > 1.0, 
or 


I < I 


r; » 


or  both. 

The  pressure  for  hot,  expanded  states  is 


Pp  = alp  + 


il£ 


I 


1 


+ Ape 


(3.2) 


where 


n < 1 .0 


nnd 


Is'- 


The  pressure  for  the  transition  state  is 


^T 


^ '’C 


(3.3) 


where 


n < 1.0 


and 


In  the  above  equations, 


I <!<!'. 
-S  s 


(3.4) 


and 


M = n - 1 . 


(3.5) 
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In  thn  above  equations,  is  the  referenee  specific  internal  energy. 

,,  IS  the  specific  internal  energy  of  the  material  as  it  is  brought  to  the 
vaporization  temperature,  I^'  is  equal  to  I^,  plus  the  latent  heat  of  vapor- 
ization for  that  material,  and  represents  the  reference  density  The 

not  represent  physical  quantit  iM . “atorial , and  do 

The  material  assumed  to  represent  the  soil  is  dry  tuff^a  with  s 

“ pr™”Hert-:r‘?h:  7e  ^ -'‘s'.!  o™ 

there'’;?rn"r:„:;\'f:  l“et:,r;:n:  fL"ty“''',?a?e'““d"“"'"’  . 

of  that  sSoil.  The  soil  qnH  i « ^ type,  state,  and  propertje 

estimate,  based  on  the  informatiorprondedl°'‘''a,;f ol!’’!'  "fr!®"*  “ '>“* 

W F TmKm-rse.2  1 a,u  . ptoviaeu  , and  on  consultations  with 

aisIime*l;St‘e™ia|vf’'|r:r'=‘"‘‘'°'  I’"-’  -nstants  for  UUs 


es 

t 


and 


a = 0..S, 

b = 1.3, 

A = 4..';  (ll’a, 

B = .3.0  GPa, 

=6.0  MJ/kg, 

a = S . 0 , 

e = s.o. 

= .3.2  MJ/kg, 

Ij."  -■  18.0  .M.J/kg, 
=1.7  Mg/m^, 


Allen,  R.  T., 
General  Atomic 

^ ^Johnson,  W.  E, 
oonmunication 


Eqmtion  of  State  of  Rooks  and  Minerals 
, Divtsion  of  General  Dijnamias,  San  Diego 


, Computer  Code  Consultants, 
(November  1974). 


Salona  Beach 


GAMD  7834, 

CA  (March  1967). 

CA,  private 
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IV.  CALCULATION  OF  TARGET  PLATE  LOADS 


As  in  the  previous  computations^ ^ the  hydrocode  computation  is 
begun  with  the  explosive  being  represented  by  a high-pressure,  quiescent 
gas  occupying  the  original  volume  of  the  explosive.  A complete  conver- 
sion of  chemical  to  internal  energy  is  assumed.  At  time  zero,  the  pres- 
sure in  the  explosive  products  as  computed  by  using  Equation  (1)  with 
Equation  (2.2)  is  11.76  GPa,  the  pressure  in  the  air  is  100  kPa,  and  the 
soil  is  at  zero  stress.  It  would  have  been  more  precise  to  begin  the 
problem  with  the  soil  at  a compressive  stress  of  100  kPa;  however,  the 
error  introduced  is  negligible:  100  kPa  is  insignificant  as  compared 
with  the  soil  stresses  caused  by  the  explosion.  Figure  2a  shows  the  pres 
sure  in  the  computational  flow  field,  plotted  against  both  the  radial  and 
axial  coordinates  at  1.0  gs.  The  plotting  program  that  produces  these 
plots  of  pressure  in  the  flow  field  plots  the  cell  pressure  at  the  upper 
right  corner  of  the  cell,  rather  than  at  the  cell  center,  the  actual  loca 
tion  of  the  pressure  value.  This  causes  the  plot  to  show  an  apparent  gap 
along  each  axis.  These  two  gaps  are  actually  the  first  column,  and  first 
row,  of  cells.  I'igure  2b  shows  a spatial  cut  of  pressure  up  the  first 
column  of  flow  field  cells,  also  at  1.0  gs.  Shocks  have  not  yet  formed 
at  this  early  time,  although  the  axial  movement  of  the  expansion  waves 
into  the  high-pressure  explosive  products  is  evident.  Figures  3a  and  3b 
show  similai'  plots  at  3.2  gs.  At  this  early  time,  shocks  are  well  devel- 
oped, the  expanding  explosive  products  initially  causing  planar  shock 
waves  to  develop  axially,  and  a cylindrical  shock  wave  to  develop  radi- 
ally. Pressure  plots  for  the  computational  flow  field  show  that  these 
shock  waves  gradually  merge  to  form  a nearly  spherical  shock  in  the  soil. 
The  shock  front  expands  unimpeded  into  the  soil  until  the  contact  surface 
between  the  soil  and  the  air  is  reached.  Analysis  of  the  DORF  computa- 
tion indicates  that  the  shock  wave  first  meets  the  soil-air  contact  sur- 
face at  17.1  gs  in  the  innermost  radial  cell,  the  point  of  intersection 
of  the  shock  wave  and  contact  surface  moving  out  radially  thereafter. 

The  incident  shock  in  this  innermost  cell  is  still  strong,  the  peak 
pressure  at  the  time  of  arrival  at  the  soil-air  contact  surface  being  5.2 
GPa.  A strong  expansion  wave  is  reflected  back  into  the  soil,  and  a 
relatively  weak  shock  wave,  whose  peak  pressure  is  on  the  order  of  1.0 
MPa,  is  transmitted  into  the  air.  As  the  point  of  intersection  of  the 
incident  shock  wave  with  the  soil-air  contact  surface  moves  radially  out- 
ward, the  soil  bulges  up  behind  it,  traveling  at  relatively  high  speed 
toward  the  target  plate. 

The  computation  indicates  that  the  computational  air  shock  first  • 
strikes  the  top  boundary  at  85  gs,  the  reflected  air  shock  in  the  inner- 
most radial  cell  at  the  top  boundary  reaching  a peak  pressure  of  4.4  MPa 
at  93  gs,  as  shown  in  Figure  4a;  the  soil  first  arrives  at  this  cell  at 
112  gs,  reaching  a peak  pressure  of  1.0  GPa  at  119  gs.  Figure  4b  shows 
the  pressure  in  the  computational  flow  field  at  128  gs,  shortly  after 
the  peak  pressure  has  been  reached  in  the  inner -most  radial  cell  at  the 
top  boundary.  The  point  of  peak  pressure  on  the  top  boundary  is  moving 
out  radially,  and  the  pressure  in  the  cells  interior  to  that  peak  is 
relieving  as  the  soil  reflects  off  the  top  boundary.  This  may  be  seen  in 
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PRESSURE  (GPo) 


AXIAL  DISTANCE  (cm) 


Figure  2b:  Spati.il  cut  of  pressure,  first 
column  of  cells  along  the  axis  of  symmetry, 
1.0  us.  No  shocks  are  evident. 
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PRESSURE  (GPo) 


Figure  3b:  Spatial  cut  of  pressure,  first  column 
of  cells  along  the  axis  of  symmetry, 

3.2  us.  Shocks  have  formed. 


Figure  4a:  The  pressure-time  history  in  the 
most  radial  cell  on  the  top  boundary. 


inner- 
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Figure  4c,  a plot  of  pressure  versus  radius  at  the  top  boundary.  In 
part,  the  jaggedness  of  this  pressure  plot  may  be  due  to  computational, 
and  not  real,  phenomena,  as  is  common  in  explicit-time-step,  Eulerian 
hydrocodes  in  which  artificial  viscosity  is  not  invoked,  as  in  the  case 
for  this  computation.  Figure  4d  shows  a spatial  cut  of  pressure  along 
the  first  column  of  cells  along  the  axis.  The  shock  wave  in  the  soil 
has  already  passed  through  the  bottom  boundary  of  the  grid.  Although  the 
bottom  boundary  is  designed  to  be  transmissive,  a shock  is  reflected 
back  up  the  grid.  The  pressure  behind  i ne  incident  shock  at  the  bottom 

boundary  is  on  the  order  of  0.6  GPa . The  reflected  wave,  at  128  ps,  is 

centered  at  an  axial  distance  of  6 cm,  with  a peak  pressure  of  0.17  GPa, 
and  is  about  to  strike  a recompression  shock  travelling  down  the  grid. 
Figure  4e  shows  a velocity  vector  map  of  the  flow  field.  The  vector 
lengths  are  scaled  relative  to  the  largest  velocity  magnitude  in  the 

flow  field  at  that  time.  Any  velocity  vector  whose  magnitude  is  less 

than  ten  per  cent  of  the  largest  magnitude  is  not  plotted.  Cratering 
caused  by  the  expanding  explosive  products  is  well  under  way.  The  low- 
velocity  toroidal  center  of  the  recompression  region,  centered  at  a 
radius  of  2.8  cm,  and  a height  of  20  cm,  is  also  evident  in  Figure  4e. 

The  hydrocode  solution  reveals  an  interesting  phenomenon  associated 
with  the  development  of  the  crater  in  the  soil.  The  soil  that  was 
originally  above  the  charge  moves  upward  until  it  encounters  the  reflect- 
ing top  boundary.  It  tnen  rebounds  back  into  the  developing  crater,  but 
its  travel  is  restricted  by  the  high  pressure  explosive  products.  Mean- 
while, the  resistance  of  the  surrounding  soil  restricts  the  radial  expan- 
sion of  the  crater.  These  axial  and  radial  constraints  lead  to  the  for- 
mation of  an  annular  jet  of  soil,  which  apparently  becomes  the  mechanism 
for  ejecting  the  soil  from  the  expanding  crater.  Figures  5a  through  5e 
show  a time-sequence  of  the  annular  soil  jet  development  and  corresponding 
crater  growth.  Those  regions  that  show  overlapping  of  contact  surfaces 
are  regions  where  mixing  of  soil,  air,  and  explosive  products  occurs. 

The  plots  of  the  contact  surfaces  are  constructed  by  connecting  tracer 
particles  that  were  placed  along  the  original  contact  surfaces  at  time 
zero.  While  the  tracer  particles  arc  fairly  reliable  in  following  the 
contact  surfaces  in  the  flow  field  regions  away  from  the  top  boundary, 
there  does  seem  to  be  some  loss  of  definition  of  the  contact  surfaces 
near  the  top  reflecting  boundary,  most  likely  due  to  the  reflections 
off  that  top  boundary.  Figures  5a  through  5e  show,  and  analysis  of  the 
computation  confirms,  that  there  are  large,  continuous  areas  of  soil 
contact  on  the  top  boundary,  and  hence  the  simulated  target  plate,  at 
all  times  after  the  initial  soil  contact  at  112  iis.  Figures  5a  through 
5e  also  show  the  development  of  the  reflected  wave  from  the  top  boundary 
that  travels  back  into  the  growing  cr.'iter.  The  impact  of  the  annular 
soil  jet  provides  the  bulk  of  the  target  plate  loading  after  130  ps.  By 
500  us  the  soil  jet's  radial  expansion  has  slowed  considerably.  During 
most  of  the  time  after  300  ps,  the  soil  jet  impacts  the  top  boundary  at 
a radius  corresponding  to  the  outer  corners  of  the  target  plate.  By 
603  ps  the  soil  jet  has  moved  off  the  area  covered  by  the  25.4  cm  square 
plate,  and  the  computation  is  halted  at  this  time.  At  this  time,  it 
appears  that  the  soil  throughout  the  lower  sections  of  the  flow  field  has 
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Figure  4c:  The  pressure  along  the 
top  boundary  at  128  ps. 
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Figure  4d:  Spatial  cut  of  pressure,  first  column  of 
cells  along  the  axis  of  symmetry,  128  ys. 
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Figure  Sa:  The  computational  flow  field  at  205  ys;  the  soil 
is  beginning  to  move  out  radially;  an  annular  soil  jet  is 
beginning  to  develop.  Velocity  vectors  are  plotted  for 
every  second  cell  in  the  axial  and  radial  directions. 
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I igure  5b:  The  computational  flow  field  at  301  ps,  reflected  wave 
from  top  boundary  developing,  annular  soil  jet  development 
continuing.  Velocity  vectors  are  plotted  for  every 
second  cell  in  the  axial  and  radial  directions. 
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Figure  5c:  The  computational  flow  field  at  398  ps,  the  reflected 
wave  from  the  top  boundary  and  the  annular  soil  jet  fully 
developed.  Velocity  vectors  are  plotted  for  every 
second  cell  in  the  axial  and  radial  directions. 
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Figure  5e:  The  computational  flow  field  at  60.'?  ps,  the  annular 
soil  jot  persisting  near  16  cm  radius,  the  reflected  wave  from 
the  top  boundary  still  moving  back  into  the  crater. 
Velocity  vectors  are  plotted  for  every  second  cell 
in  the  axial  and  radial  directions. 
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begun  a bulk  upward  movement.  This  is  almost  surely  an  artifact  of  the 
code,  as  will  be  discussed  later.  Figure  6 shows  the  pressure  distribu- 
tion on  the  top  boundary  of  the  computational  flow  field  at  398  ys.  The 
large  pressure  peak  is  due  to  the  soil  jet  (see  Figure  5c).  This  type 
of  pressure  distribution  is  typical  for  the  time  during  which  the  soil 
jet  exists. 

As  mentioned  previously,  the  computation  is  carried  out  to  603  ps, 
real-time.  The  target  plate  is  represented  by  a rigid,  stationary, 
square  plate,  25.4  cm  on  a side,  centered  above  the  mine.  The  total  im- 
pulse, gage,  on  the  target  plate  is  computed  by: 


Total  Impulse  = 


where 


J = the  number  of  points  in  time  at  which  loading  data  for 

the  top  boundary  was  stored  on  a data  tape  (loading  data 
was  stored  every  5 computational  cycles  for  this  problem). 


N 

A. 


P.  . 
i.J 


the  number  of  flow  field  cells  intersected  by  the  target 
plate , 

the  area  of  the  intersection  of  the  i^^  flow  field  cell 
with  the  target  plate  (please  see  following  discussion 
for  details) , 

the  pressure  in  the  i'’*^  flow  field  cell  at  the  j^^'  point 
in  time, 


P = the  assumed,  constant,  atmospheric  pressure  behind  the 
rigid,  immobile  target  plate,  and 

t = the  time. 


Figure  7 shows  a top  view  of  the  grid  for  the  hydrodynamic  flow  field 
used  in  this  DORF  computation,  with  the  axis  of  symmetry  normal  to  the 
page  and  located  at  point  0.  The  square  ABCD  represents  the  square  tar- 
get plate,  where  each  side  is  25.4  cm.  The  concentric  circles  indicate 
the  outer  radial  boundary  of  every  fifth  flow  field  cell  in  the  radial, 
or  I,  direction.  The  grid  spacing  for  the  deformation  computation  (see 
Section  V),  20  even  cells  per  side  of  square  AEOF,  is  also  shown.  The 
loading  as  computed  by  using  Equation  (5)  includes  that  contribution 
from  each  flow  field  cell  which  intersects  the  square  ABCD,  proportional 
by  the  intersect  of  the  flow  field  cell  with  the  target  plate. 
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PRESSURE  (MPo) 


^Igure  6.  Pressure  distribution  on  top  boundary  of 
computational  flow  field  at  .^98  ps  11"^^ 
pressure  peak  due  to  soil  jet  impact. 
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Figure  7.  A top  view  of  the  intersection  of  the  25.4  cm 
square  target  plate  ABCD  with  the  finite  difference 
grid  for  the  DORF  loading  computation,  and  the 
grid  for  the  RF.PSIL  response  computation 
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As  mentioned  previously,  the  target  is  represented  by  a rigid  sta- 
tionary, square  plate,  25.4  cm  on  a side,  centered  above  the  mine. 

Figure  8a  shows  a plot  of  total  impulse,  gage,  versus  time  on  the  target 
plate.  As  may  be  seen  in  Figure  8a,  the  initial  airblast  makes  an  insig- 
nificant contribution  to  the  total  impulse.  The  subsequent  total  impulse 
is  almost  entirely  due  to  the  impact  of  the  soil.  The  total  impulse 
reaches  an  apparent  plateau  at  0.77  kN-s  at  603  ps,  4,5"o  above  that  pre- 
dicted by  the  correlation  function  generated  by  Westine  and  Hokanson^*’. 
The  impulse  correlation  function  by  Westine  and  llokanson^®  provides  only 
total  values  for  impulse,  and  not  a time-history.  Figure  8b  is  a plot 
of  resultant  force,  gage,  versus  time  on  the  target  plate,  computed  by 


i = N 

Force  = ^ A.  (P.  - 1'.^) , 
i-1 


(6) 


where 


P. 

1 


is  the  pressure  in  the  i 


th 


flow  field  cell, 


and  the  other  variables  are  as  defined  for  E<iuation  (5).  The  large, 
initial  peak  corresponds  to  the  arrival  of  the  soil  at  the  target  plate. 


The  DORF  computer  simulation  also  furnishes  the  time-history  of 
the  crater  volume,  as  may  be  seen  in  Figure  9.  The  crater  undergoes 
rapid,  monotonic  growth  up  until  325  ps,  at  which  time  the  growth  becomes 
somewhat  erratic.  The  peak  crater  volume  of  4.52  xlO^  cm^  occurs  at  400 
ps;  by  603  ps  it  has  fallen  to  4.02  xlO^  cm^  This  decrease  in  volume 
at  later  time  is  most  likely  an  artifact  of  the  code,  possibly  caused  by 
the  artificial  reflection  of  the  shock  wave  back  into  the  soil  at  the 
bottom  transmissive  boundary.  It  appears  that  sometliing,  most  likely 
the  reflected  wave  mentioned  above,  is  causing  a bulk,  upward  movement 
of  the  soil.  Figures  5a  through  5e  show  the  devclojvncnt  of  this  upward 
movement.  While  this  casts  some  doubt  on  the  calculation  of  the  late- 
time crater  volume,  it  does  not  appear  to  have  affected  the  target  plate 
loading  calculation.  This  is  so  because  this  artificial  wave  was  detected 
during  one  of  the  several  routine,  intermediate  analyses  of  the  solution 
at  a time  sliortly  after  its  appearance  at  the  bottom  boundary  as  an  ap- 
parent artificial  reflection  of  the  incident  shock  wave.  The  computation 
was  stopped  at  603  ps,  when  the  leading  edge  of  this  wave  was  still  sev- 
eral cells  away  from  the  top  boundary,  the  area  of  primary  interest  in 
this  computation. 
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figure  8a:  Impulse,  gage,  integrated  over  space  and 
time  on  the  25.4  cm  square  target  plate. 
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In  a previous  publication  by  Westine^^,  and  later  in  reports  by 
Westine,  Baker,  and  Wenzel^'*  and  Wenzel  and  bsparza^"*,  the  volume  for 
the  crater  made  by  a mine  with  no  target  above  is  predicted  by  the  rela- 
tion 


for 


and 


V = 8.49 


id) 


0.66 


>0.4 


(7) 


d > 0.0 


Here,  V represents  volume  in  cubic  feet,  W represents  charge  weight  in 
pounds,  and  d represents  the  aepth  of  burial,  in  feet  from  the  overburden 
surface  to  the  top  of  the  charge.  Westine^^  indicates  that  the  crater 
volume  is  relatively  independent  of  the  overburden  material.  For  the 
case  reported  here,  d is  0.188  feet  and  W is  0..883  Ibm  (therefore 
7 /24 

W /d  = 4,.S4  > 0.4]  and  F.quation  (7)  predicts  a crater  volume  of  1,95 
cubic  feet,  or  5.5.5  cm^,  12. J times  greater  than  that  predicted  by 

the  DORF  computation  reported  here.  The  reason  for  this  large  discrepancy 
is  not  known,  although  it  may  be  that  liquation  (7)  is  not  applicable  for 
small  charges.  The  smallest  charge  weight  apparent Iv  used  to  generate 
the  correlation  function  is  116  kg  of  TNi'.  Westine'^  also  jiresents  a corre- 
lation for  crater  radius  in  the  form  of  a plot,  which,  for  the  case  pre- 
sented here,  predicts  a crater  radius  ratio  of  R/d  = 8. .52.  Because  d is 
5.74  cm,  the  predicted  crater  radius,  assumed  to  be  measured  at  the  orig- 
inal ground  level,  is  17.7  cm.  The  correspoiul  i ng  crater  radius  predicted 
by  the  DORF  calculation  is  1.5.5  cm,  so  the  cr.itio’  r.idius  [redicted  using 
Westine's  correlation  is  .5.5  times  that  predicted  lyv  the  IKIRF  calculation, 
which  is  at  least  consistent  with  the  volume  il  i screp.ancy . Aiso,  the 
crater  radius  and  volume  comi)uted  from  the  PORI  solution  .are  for  the  time 
603  p.s,  the  time  at  which  the  computation  was  stopped.  At  this  time,  the 
soil  velocity  in  the  annular  soil  jet  was  on  the  order  of  230  m/s.  It  is 


r.  i'.  , crinj , > •_/’  Tt  I'Vivriticliuni  cc, , 

Penjiirnrnon  PreUfi,  l.'>70,  17’/,  2,  No.  .V’. 
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assumed  that  crater  growth  would  have  continued  beyond  this  time  if  it 
had  not  been  necessary  to  stop  the  computation  because  of  the  strong 
artificial  wave  coming  from  the  bottom  boundary. 


V.  TARGET  PLATE  RESPONSE 

The  REPSIL^’  (Response  of  Elastic-Plastic  Shells  to  Impulsive  Loads) 
structural  response  program,  selected  to  predict  plate  response  to  the 
buried  land  mine  explosion,  is  capable  of  treating  large  transient  defor- 
mations. The  program  accepts  the  loads  generated  by  DORF  on  the  target 
plate  as  a function  of  time  and  spatial  coordinate.  These  features  make 
REPSIL  suitable  for  predicting  plate  response  to  a buried  land  mine  explo- 
sion . 


The  computational  model  of  the  target  plate  takes  advantage  of  the 
dual  symmetry  of  the  problem;  only  the  lower  left  quarter  of  the  target 
plate,  as  shown  in  Figure  10,  is  treated.  Although  the  SwRI  experimental 
boundary  conditions  did  not  achieve  an  ideally  clamped  edge  condition, 
modeling  them  as  such  in  REPSIL  was  the  best  available  approximation  to 
the  experimental  conditions.  A 20  x 20  finite  difference  mesh  is  used  in 
modeling  one-quarter  of  the  target  plate. 

The  target  plate  material  is  reported'^  to  be  armor  plate  with  a 
density  of  7.833  Mg/m^  and  a yield  strength  of  789.45  MPa.  Rolled  homo- 
geneous armor  (RHA)  was  assumed  for  other  material  properties  needed  by 
REPSIL  and  are  listed  below; 

Young's  modulus  E * 204.3  GPa, 

Poisson's  ratio  v = 0.3, 

Mass  density  p = 7.833  Mg/m^ . 

In  the  plastic  range  negligible  error  is  introduced  by  assuming  that  RhA 
strain  hardens  independent  of  the  strain  rate^®.  This  behavior  is  modeled 
in  REPSIL  by  using  a polygonal  approximation  to  the  uniaxial  loading 
curve.  The  slope  of  the  polygonal  approximation  for  RHA  changes  at  the 
following  stress-strain  points: 

= 789.45  MPa  Cj  = 0.004 

= 919.76  MPa  = 0.025 

Oj  = 962.58  MPa  = 0.075. 


^^Benokt  R.  F.  and  RohitailXe,  J.  L. , "Teneile  Sti'sea-Strain  Curves — 111, 
Foiled  Homogeneous  Armor  at  a Strain  Rate  of  0.42  3~^,"  BRL-MR-2760,  US 
Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD  (June  1977), 
(AD  #A041560) 
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The  dynamic  loading  of  the  rigid  target  plate  determined  by  the  DORF 
code  provides  the  loading  function  P(r,  t)  for  the  structural  response 
computations.  These  pressures  are  linearly  interpolated  in  time  and  I'adius 
from  the  charge  center  to  obtain  the  load  at  a given  mesh  point  for  each 
cycle  of  calculation  in  REPSIL  (see  Figure  7). 

The  deformation  patteni  of  the  target  plate  predicted  by  REPSIL  is 
symmetric  and  free  of  any  local  denting,  as  can  be  seen  in  Figures  11  and 
12,  which  show  the  target  plate  at  275  ys  and  575  Ms  after  detonation 
of  the  land  mine.  Local  denting  is  not  predicted  due  to  the  numerical 
approximation  of  the  loading;  pressure  spikes  narrower  than  one  mesh 
spacing  are  not  treated. 

The  jagged  plate  contour  near  the  clamped  euge  in  Figure  12  is  due 
to  the  numerical  modeling  of  a clamped  edge.  Near  the  clamped  edge  a 
one-sided  differencing  scheme  is  used  to  approximate  derivatives.  At  the 
clamped  boundary  the  components  of  the  normal  must  remain  fixed  at  their 
original  values  and  hence  the  slope  is  also  fixed  at  its  original  value 
of  zero.  Due  to  the  one-sided  difference  formula  for  the  first  deriva- 
tive used  in  REPS  1 1., 


dll 
n 

d\ 


5 u"  - 4 u"  , . u'\  , 
m m+1  m+2 


(u”  = 0 at  the  boundary) , 


the  normal  component  of  displacement  at  one  mesh  spacing  in  fri'm  the 
clamped  edge,  u'’  , , is  1/4  the  value  of  the  corresponding  normal  cor.ipo- 
ni+l  . , n 

nent  of  displacement  two  incsii  r.oacings  in  from  the  clamped  edge, 

This,  coupled  with  the  central  differencing  to  obtain  first  order 
derivatives  away  from  the  clamped  edge,  leads  to  the  jagged  euge  condi- 
tion near  the  clamped  edge. 

The  initial  soil  impact  with  the  target  plate  is  indicated  on 
Figure  15,  which  shows  the  total  deflection  history  for  the  midpoint  of 
the  target  plate.  As  anticipated,  the  maximum  residual  deformation  occurs 
at  the  midpoint  of  the  target  plate.  REPSIL  calculates  a maximum  residual 
deflection  of  5.59  cm  when  the  target  plate  reaches  its  final  equilibrium 
configuration  (kinetic  energy  of  plate  is  zero)  as  shown  in  Figure  13. 

The  amount  of  energy  delivered  to  the  target  plate  from  a buried 
land  mine  explosion,  based  on  the  rigid  body  motion  of  a mechanical 
system  having  a single  degree-of- freedom' ® , is  the  square  of  the  total 
impulse  delivered  by  the  explosive  products  and  soil,  divided  by  twice 
the  plate  mass:  this  compares  well  with  the  total  energy  in  the  system 
as  computed  by  REPSIL.  The  total  impulse  delivered  to  the  target  plate 
as  computed  using  the  DORF  hydrocode  is  0.77  kN-s,  yielding  an  energy 
of  86.47  kJ  being  delivered  to  the  target  plate.  The  total  impulse 
delivered  to  the  target  plate  using  the  empirical  relationship  developed 
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Figure  11,  Target  plate  deflection,  time  = 275  us 
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Figure  1.^.  Deflection  time-history  of  the 


by  Westine  and  Hokanson*®  is  0.741  kN-s,  yielding  an  energy  of  80.43  kJ. 
The  amount  of  energy  delivered  to  the  target  plate,  as  computed  by  REPSIL, 
is  the  external  work  done  on  each  mesh  integrated  over  the  plate  and  is 
represented  by  the  external  work  curve  (see  Figure  14)  which  represents 
the  work  done  on  the  target  plate  by  the  external  loads.  The  88.41  kJ 
computed  by  REPSIL  compares  well  with  the  energy  values  computed  above. 

As  can  be  seen  in  Figure  14,  most  of  the  energy  delivered  to  the  target 
plate  is  dissipated  in  plastic  deformation  of  the  target  plate,  resulting 
in  a deformation  profile  that  agrees  well  with  predictions  from  empirical 
relationships.  In  the  SwRI  experiments*®,  a portion  of  the  energy  de- 
livered to  the  target  plate  is  dissipated  in  accelerating  the  target  plate 
and  collar  assembly. 

Only  the  maximum  residual  deformation  of  the  target  plate  is  reported 
in  Reference  16  for  test  No.  36.  However,  a more  detailed  comparison  of 
the  REPSIL  predictions  can  be  made  by  using  the  empirical  expression 
developed  by  Westine  and  Hokanson*®  for  predicting  the  plate  deformation 
profile  for  a clamped- clamped  flat  plate  centered  over  a buried  land  mine. 
The  expression  provides  a cubic  fit  to  existing,  experimentally  observed, 
deflection  profiles  for  flat  target  plates  subjected  to  a centered  mine 
blast  and  is  restated  below: 


w 


[■ 


(9) 


where, 

Wq  is  the  maximum  plate  deflection, 

w is  the  observed  deflection  at  point  (x,  y) , 

X is  the  plate  coordinate  in  X direction,  measured  from  the 
charge  center, 

y is  the  plate  coordinate  in  Y direction,  measured  from  the 
charge  center, 

X is  the  half  span  of  plate  in  X direction,  and 

Y is  the  half  span  of  plate  in  Y direction. 

Using  Equation  (9),  a predicted  deformation  profile  along  either  symmetry 
plane  of  the  quarter  plate  has  been  determined  and  is  plotted  in  Figure 
15;  REPSIL's  predicted  plate  deformation  profile  is  also  plotted.  Both 
curves  are  scaled  to  the  maximum  residual  plate  deformation  of  5.131  cm 
which  was  predicted  by  the  Westine  and  Hokanson*®  procedure.  The  excel- 
lent correlation  between  the  curves  illustrates  that  the  Westine  and 
Hokanson  procedure  provides  a quick  estimate  for  the  maximum  residual 
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structural  response  of  target  plate. 


scaled  deflection 


*/X,  SCALED  DISTANCE  FROM  PLATE  CENTER 


Figure  IS.  Scaled  target  plate  deformation 
profile  along  either  symmetry  plane. 


deformation  and  plate  profile.  However,  by  using  the  loading  history  for 
the  target  plate  generated  by  the  WR}’  code  as  the  loading  function,  the 
RiiRSlL  code  provides  a detailed  analysis  of  the  target  plate,  including 
such  items  as  both  elastic  and  plastic  strain  predictions  at  predeter- 
mined locations,  target  plate  strain  energy,  and  graphical  representation 
of  the  response. 

From  Figures  8a  and  11  it  can  he  seen  that  a portion  of  the  impulse 
is  delivered  to  the  target  plate  after  the  external  work  has  leveled  off. 
This  may  cast  some  doubt  on  the  assumption  that  the  OORF  calculations  can 
be  uncoupled  from  the  RFPSIL  calculations.  This  can  be  resolved  by  ob- 
serving in  Figure  8b  that  the  resultant  force  on  the  target  plate  decays 
rather  rapidly,  which  allows  for  little  plate  deformation  to  occur  even 
though  contributions  arc  being  made  to  the  iminilse. 


VI.  SUMMARY 

The  numerical  simulation,  reported  here,  of  a field  experiment^® 
performed  and  reported  by  SwRI  has  been  compared  with  empirical  correla- 
tion functions'®.  The  DORF  hydrocode  solution  for  the  total  impulse 
delivered  to  the  rolled  homogeneous  armor  target  plate  is  4.5%  above 
that  computed  using  the  SwRI  correlation  function  for  impulse.  The  maxi- 
mum residual  deformation  profile  of  the  target  plate,  predicted  by  the 
REPSIL  structural  response  program,  agrees  to  within  10%  of  the  profile 
computed  using  the  SwR!  correlation  function.  The  DORF  hydrocode  and  the 
REPSIE  structural  response  code  give  a reasonable  estimate  for  the  plate 
deflection,  even  though  the  loading  and  response  calculations  arc  un- 
coupled. 


VI  I.  RECOMMENDATIONS 

Although  the  above  stated  comparisons  of  computed  results  versus 
the  empirical  correlation  functions  show  excellent  agreement,  there  is 
considerable  room  for  improvement  in  both  the  DORF  code  and  in  the  tech- 
niques used  here.  The  more  significant  recommendations  for  improvement 
and  for  additional  efforts  are  as  listed  below. 

A.  Add  an  accurate  explosive  detonation  routine  to  the  version  of 
DORF  used  here. 

B.  Add  an  appropriate  equation  of  state  to  DORF  to  model  the  explo- 
sive products. 

C.  Improve  the  DORF  code's  ability  to  model  soil  and  investigate 
the  apparent  crater  volume  and  radius  discrepancy. 

D.  Improve  the  numerical  formulation  of  the  transmissive  boundaries 
in  DORF,  with  particular  attention  to  the  formulation  of  the  bottom  boun- 
dary . 


E.  Modify  the  differencing  technique  In  DORP  to  eliminate  the  errors 
caused  by  DORF's  equal  weighing  of  unequal  cell  volumes. 

F.  Establish  a more  detailed  understanding  of  the  artificial  visco- 
sity option  in  DORF.  It  was  not  used  in  this  computation  because  it 
seems  to  have  a significant  effect  on  the  solution  beyond  simply  improving 
numerical  stability,  which  is  the  all  too  generally  accepted  reason  for 
its  use. 


G.  Develop  a clamped  boundary  model  for  REPSIL  which  allows  plate 
slippage. 

H.  Couple  a hydrodynamic  code  to  compute  loading  with  a structural 
response  code  suitable  for  computing  plate  response. 
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